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A coexistent phase of spin polarization and color superconductivity in high-density 
QCD is investigated using a self-consistent mean-field method at zero temperature. The 
axial- vector self-energy stemming from the Fock exchange term of the one- gluon- exchange 
interaction has a central role to cause spin polarization. The magnitude of spin polar- 
ization is determined by the coupled Schwinger-Dyson equations with a superconducting 
gap function. As a significant feature, the Fermi surface is deformed by the axial-vector 
self-energy and then rotation symmetry is spontaneously broken down. The gap function 
results in being anisotropic in the momentum space in accordance with the deformation. 
As a result of numerical calculations, it is found that spin polarization barely conflicts 
with color superconductivity, but almost coexists with it. 



I. INTRODUCTION 

Recently much interest is given for high-density QCD, especially for quark Cooper-pair condensa- 
tion phenomena at high-density quark matter (called as color superconductivity (CSC)), in connection 
with, e.g., physics of heavy ion collisions and neutron stars [1-3]. Its mechanism is similar to the BCS 
theory for the electron- phonon system [4], in which the attractive interaction of electrons is provided 
by phonon exchange and causes the Cooper instability near the Fermi surface. As for quark matter, 
the quark-quark interaction is mediated by colored gluons, and is often approximated by some effective 
interactions, e.g., the one-gluon-exchange (OGE) or the instanton-induced interaction, both of which 
give rise to the attractive quark-quark interaction in the color anti-symmetric 3* channel. CSC leads to 
spontaneous symmetry breaking of color SU{i) into SU{2) as a result of condensation of quark Cooper 
pairs [2,3]. 

In this paper we would like to address another phenomenon expected in quark matter: spin polariza- 
tion or ferromagnetism of quark matter. We examine the possibility of the spin-polarized phase with CSC 
in quark matter. As far as we know, interplay between the color superconducting phase and other phases 
characterized by the non- vanishing mean fields of the spinor bilinears (ipTip) has not been explored except 
for the case of chiral symmetry breaking [5] . Our main concern here is to investigate the possibility of 
the quark Cooper instability under the axial-vector mean-field, (^'7^75^') which is responsible for spin 
polarization of quark matter. It would be worth mentioning in this context that ferromagnetism (or spin 
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polarization) and superconductivity are fundamental concepts in condensed matter physics, and their 
coexistent phase has been discussed for a long time [6]. As a recent progress, a superconducting phase 
have been discovered in ferromagnetic materials and many efforts have been made to understand the 
coexisting mechanism [7]. 

Besides being interesting in its own right, the coexistence problem may be related to some physical 
phenomena. Recently, a new type of neutron stars, called as "magnetars", with a super strong magnetic 
field of ^ O(IO^'^G) has been discovered [8,9]. they may raise an interesting question for the origin of the 
magnetic field in compact stars, since its strength is too large to regard it as a successor from progenitor 
stars, unlike canonical neutron stars [10]. Since hadronic matter spreads over inside neutron stars beyond 
the nuclear density (po^0.16fm~'^), it should be interesting to consider the microscopic origin of the 
magnetic field in magnetars. In this context, a possibility of ferromagnetism in quark matter due to 
the OGE interaction has been suggested by one of the authors (T.T.) within a variational framework 
[11]; a competition between the kinetic and the Fock exchange energies gives rise to spin polarization, 
similarly to Bloch's idea for itinerant electrons. Salient features of spin polarization in the relativistic 
system are also discussed in Ref. [11]. Thus, it might be also interesting to examine the possibility of the 
spin-polarized phase with CSC in quark matter, in connection with magnetars. 

We investigate spin polarization in the color superconducting phase by a self-consistent framework, 
in which quark Cooper pairs are formed under the axial-vector mean-field. We shall see that this phe- 
nomenon is a manifestation of spontaneous breaking of both color SU{3) and rotation symmetries. 

We adopt here the OGE interaction as an effective quark-quark interaction. Since the Fermi mo- 
mentum is very large at high density, asymptotic freedom of QCD implies that the interaction between 
quarks is very weak [12]. So it may be reasonable to think that the OGE interaction has a dominant 
contribution for the quark-quark interaction. In the framework of relativistic mean-field theories, the 
axial- vector and tensor mean-fields, which stem from the Fock exchange terms, ('!^757^V') and (^"'/ii^V') i 
may have a central role to split the degenerate single-particle energies of the two spin states, and then 
leads to spin polarization, e.g., see [13] for discussion in nuclear matter. As for quark matter, several types 
of the color singlet mean-fields appear after the Fierz transformation in the Fock exchange terms, but 
we retain only the axial- vector mean-field as the origin of spin polarization, because the OGE interaction 
by no means holds the tensor mean-field due to chiral symmetry in QCD, unlike nuclear matter [13]. 
Presence of the axial-vector mean-field deforms the quark Fermi seas according to their spin degrees of 
freedom, and thereby the gap function should be no more isotropic in the momentum space. We assume 
here an anisotropic gap function A on the Fermi surface by a physical consideration and solve the cou- 
pled Schwinger-Dyson equations self-consistently by way of the Nambu formalism to find the axial-vector 
mean-field Ua and the superconducting gap function A. Thus we discuss the interplay between spin 
polarization and superconductivity in quark matter. 

In Section 2 we give a framework to deal with the present subject. The explicit structure of the 
anisotropic gap function A in the color, flavor, and Dirac spaces is carefully discussed there and in the 
Appendix B and Appendix C. Numerical results about Ua and A are given in Section 3, where phase 
diagram of spin polarization and color superconductivity is given in the mass-baryon number density 
plane. Section 4 is devoted to summary and concluding remarks. 
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II. FORMALISM 



In this section we present our formalism to treat CSC and spin polarization. We consider quark 
matter with flavor SU{2) and color SU (3) symmetries, and assume that the interaction action is described 
by the OGE interaction as 



(1) 



are the 



where t/' is the quark field, -D^^(a;, y) is the gauge boson (gluon) propagator, and Xa = 1,2,- 
SU{3) Gell-Mann matrices. Using the Nambu formalism [2,14] the effective action is given within the 
mean-field approximation as 



Imf = \ 




g-Hp) 




(2) 



(3) 



with the inverse quark Green function 

where /i = 70/i with the chemical potential fi. F is a self-energy and A is the gap function for the quark 
Cooper pair; both terms V and A should be provided by the Fock exchange terms of the OGE interaction. 
We define here tjjc{k) and V as 

Mk) = C^^(-fc), (4) 
V = CV^C-^ (5) 

with the charge conjugation matrix C which is explicitly given by 17270 in Dirac representation. 
The Green function G{p) can be written straightforwardly from eq.(3) as 

Gn(p) Gi2(p)\ 

G2l{p) G22(P) / 



G{P) 



with 



Glib) =[^-m + ^ + V{p) - 7oA(p)^7o {i)-m-fi + V{p)) ^ A{p) 
G2M = -{i>-m- ^ A(p)Gii(p). 



(6) 

(7) 
(8) 



Following Nambu's argument [14], we impose the self-consistency condition to obtain the Hartree-Fock 
ground state such that the self-energy by the residual interaction, Sijes. , vanishes: 



Sfies. = ^M.F. — ^Int. — 0, 



where Em.f. is defined by 

^M.F.{k)=Go\k)-G-\k) 
{i>-m + /i)-^ 



with G(y{p) = 



V{k) 7oAt(fc)7o 
A(A;) V{k) 









{]/) — m — i/i) ^ 



(9) 

(10) 
(11) 



3 



and Tijnt. is given by the use of the OGE interaction. Within the first-order approximation in g'^, S/„t. 
renders 



^int.{k)=g'J Al_D'^'>^k-p)raG{p)h (12) 

(13) 



\ _ / 

C(7''^)'^C-M \ -7^^ 



which is nothing else but the Fock exchange energy by the OGE interaction. Using eqs. (9) - (12), 
we obtain the self-consistent equation for V{k) by the use of the diagonal component of the full Green 
function (7): 

-V{k) = {-igf j -^[-ii5M^(fc-p)h^^[-iGn(p)]7.^- (14) 

The gap equation is also obtained from the off-diagonal component as 

-A(/c) = {-igf j ^Hi)M.(fc-p)]7^^(^[-iG2i(p)]7.^. (15) 

In the following sections, we present explicit forms of V{p) and A(p) and then solve their coupled equations 
(14) and (15). 

A. Fermion propagator under the axial-vector self-energy 



We, hereafter, take the static approximation for the gauge-boson propagator as 

« (16) 

where M is an effective gauge boson mass originated from the Debye screening ~ Nfg^ij,'^/{2n^) [15]. 

Since typical momentum transfer |q| at high density is of the order of the chemical potential, we may 
further introduce the zero-range approximation [16] for the propagator as 

with a typical momentum scale Q of 0(/i). This approximation corresponds to the Stoner model [17], 
which is popular in solid-state physics, and stands on the same concept of the NJL model [18] as well. 

To proceed, we assume, without loss of generality, that total spin expectation value is oriented to 
the negative ^;-direction in the spin-polarized phase which is caused by the finite axial-vector mean-field 
along the 2;-axis ^ . As shown in Ref. [13], rotation symmetry is spontaneously broken down in this phase 
while axial symmetry around the ^-axis is preserved. Then two Fermi seas of the different spin states are 
deformed accordingly. 



'^We shall see that only the space component of the axial-vector mean-field is responsible for spin polarization. 
We,hereafter, take its direction along the z-axia without loss of generality. 
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Applying the Fierz transformation for the Fock exchange energy term (14) we can see that there appear 
the color-singlet scalar, pseudoscalar, vector and axial- vector self-energies (Appendix D). In general we 
must take into account these self-energies in V, V = Us + i'^hUps + liJiU^ + l^lbU^^ with the mean-fields 
Ua- Here we introduce an ansatz: the Femri distribution holds the reflection symmetry with respect to 
the Px—Py plane, and only the mean-field parts Us, U^ and f/„^ are retained in V. Later we will see that 
the self-consistent solution is obtained with the zero-range approximation (17) under this ansatz. 

In this paper, furthermore, we disregard the scalar mean-field Us and the time component of the 
vector mean-field U^ for simplicity since they are irrelevant for the spin degree of freedom; U^ has only a 
role to shift the total energy ot the chemical potential, and may not affect any other physical properties. 
On the other hand. Us may significantly influence the spin-polarization properties through changing the 
quark effective mass. Instead of introducing the scalar mean-field explicitly, however, we treat the quark 
mass as a variable parameter, and discuss its effect in the next section. 

According to the above assumptions and considerations the self-energy V in eq.(3) renders 

1^ = 7375«7a, Ua = UI, (18) 

with the axial- vector mean-field Ua- Then the diagonal component of the Green function Gii(p) is 
written as 

Gii(p) = [g-/ - 7oAt7oG^A] (19) 

with 

G-A\p)=i>-m + ii-^^-f^UA, (20) 
G-,\p)=i>-m-ii-WllUA, (21) 

where 7573 = 7573 and Ga{p) is the Green function with the axial- vector mean- field Ua which is deter- 
mined self-consistently by way of eq. (14). 

Before constructing the gap function A, we first find the single-particle spectra and their eigenspinors 
in the absence of A, which is achieved by diagonalization of the operator G^. In the usual case of no 
spin polarization this procedure gives nothing but the free energy spectra and plane waves. Then we 
choose a gap structure on the basis of a physical consideration as in the usual BCS theory. 

From the condition that det G^^(po)U=o = one can obtain four single-particle energies e± (positive 
energies) and — e± (negative energies), which are given as 



(p) = ^f^+u\ + ^^^±2^u^/^f^,, (22) 



where the sign factor ±1 being in front of Ua indicates the energy splitting between different spin states 
due to the presence of the axial-vector self-energy, which corresponds to the exchange splitting in the 
non-relativistic electron system [17]. In the following, we call the "spin" -up (-down) states for the states 
±e+ (±6-). Eq. (22) also shows that each Fermi sea for the "spin" -up (-down) state should undergo a 
deformation and lose rotation symmetry, once Ua is finite. This is a genuine relativistic effect [13]; actually 
the exchange splitting never produces deformation of the Fermi sea in the non-relativistic ferromagnetism, 
e.g. in the Stoner model [17]. 
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Here, it would be interesting to see the peculiarities of the quark Fermi seas in the presence of the 
axial- vector self-energy. In Fig. 1 we sketch the profile of the Fermi seas projected onto the pz ~pt plane 
(Pt — \Jpi foi' the cases of (a) Ua < m, (b) Ua > m and (c) m — 0. As is already mentioned 

these seas still hold the axial symmetry around the z-axis and the reflection symmetry with respect to 
the px — Py plane. The region surrounded by the outer line show the Fermi sea of "spin" -down quarks, 
and the shaded region is that of "spin" -up quarks. 

(b) 








1 ^'""""'^-l^ 






FIG. 1. Illustrations of the Fermi surfaces in the pt — Pz plane with pt = \/p% + Py- (a) : for fi > Ua + m 
and Ua < rn. Outer closed curve corresponds to the Fermi surface of the "spin"-down state with single-particle 
energy e-(p) and inner one surrounding the shaded area to the "spin"-up state with e+(p). (b) : the same for 
(a) but Ua > m. A pair of white circles connected by dashed line represents the Cooper pair characterized by 
B„ (27). Each particle in the Cooper pair has a different color and flavor, (c) : the Fermi surfaces for m — > 0. 
The outer (inner) contour represents the Fermi surface for £_ (e+). 



We can see in Fig. la that the Fermi seas for the "spin" -down and "spin" -up states are deformed 
in the prolate and oblate shapes, respectively, where the minimum of the single-particle energy still 
resides at the origin p — 0. When Ua > to as shown in Fig. lb, there appear two minima at the points 
{pt,Pz) — (0, ±-\/?7^ — m?) for the "spin"-down quark. Hence in the massless limit, m ^ 0, the Fermi 
sea is described by two identical spheres with radii /i in the momentum space, which are centered at the 
points (pt,Pz) = {0,±Ua) (see Fig. Ic). 

In what follows we use subscript 'n'(= 1,2,3,4) for notational convenience as e„ which means 
{ei, £2, £3, €4} — {e_, e+, — e_, — £+}. We define the spinor 0„(p) that satisfies the equation G^{po = 
En — P')4'n{p) = 0, which corrcsponds to the eigenspinor with the single-particle energy e„ in the absence 
of the quark Cooper pairing. The spinor (j)n{p) is explicitly given as 



v 



(e„ - (-l)"/3p - Ua){Px - iPy)Pz 

-((-l)«/3p+TO)p2 

[-((-l)"/3p -f TO)(e„ - TO - Ua) + pI] {p. - ipy) 



\ 



pIpz 



(23) 



where Nn= ^[/5p - (-1)"to] [e„ + Ua + (-l)"/3p] /{enPp)/{2p'fpz) and Pp = ^/pl + ^. It is to be noted 
that the spinors (/>„ do not return to the eigenspinors of spin operator az even when Ua ^ 0, but become 
mixtures of them, see Appendix A. Introducing the projection operator A„ — (j>n(t>n with properties 
and J2n A„ = 1, we can recast Ga{p) in the spectral representation into 
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G^ = J2 ^^70 (24) 

= ^{po - e„ + /x)7oA„. (25) 



B. Gap structure 

In this subsection, we give the exphcit form of the gap function A in the Dirac, color and flavor 
spaces, and then calculate the diagonal component of the full Green function Gii{p) in eq. (7), provided 
that only the axial-vector self-energy is taken for V{p) in eq. (14). In general various types of the gap 
structures are possible in the Dirac, color and flavor spaces; they depend on the form of interaction and 
the quark mass [2,19], especially on the strange quark mass [20]. Here we suppose a simple gap structure 
from a physical consideration, disregarding the finite mass eff'ect. 

Using the the spinor ^„(p) we assume that the gap function A in eq. (15) has a following form in the 
color and flavor spaces: 

A(p) = ^A„(p)i?„(p) (26) 



with the operator Bn{p), 



Bnip) = 70<^-n(p)4(P)- (27) 



where the subscript '— n'(= —1,-2,— 3, —4) indicates that the single-particle energy in the spinor is 
replaced by that of opposite sign, e_„ = — e„, without change of "spin". 

One can easily see what kind of quark pairs the gap function A (26) represents. Utilizing the property, 
^'^^,{—p)C'yo(f)n{p)ocSn'n, One can find for the general spinor V(p) = ^n(p)</'n(p) with arbitrary 
coefficients a„, 

tpcBntp = ilF {-p)C^o4>-n(t>n'4^{p) oc a„(-p)a„(i?). (28) 

This equation clearly shows that two quarks included in the Cooper pairing have opposite momenta to 
each other and belong to the same energy eigenstate as illustrated in Fig. lb. 

Now we should note that the antisymmetric nature of the fermion self-energy imposes a constraint on 
the gap function [2,21], 

C7A(p)C-i = A^(-p). (29) 

Since B„ satisfies the relation CBn{p)C~^ = B'^{—p), A„(p) must be a symmetric matrix in the spaces 
of internal degrees of freedom. Taking into account the property that the most attractive channel of 
the OGE interaction is the color antisymmetric 3* one, it must be the flavor singlet state. Thus we can 
choose the form of the gap function as 

A„(p) = ea03eij^n{p), (30) 

-I a/3, ij 
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where (a/3) and {ij) are indices in three-color and two-flavor spaces, respectively. The form of gap 
function (30) in the color and flavor spaces is familiar for two-flavor CSC [2,3]. 

Using the properties of A„(p) and Bn{p), we then obtain an explicit form of Gn{p) as 

\ -1 



[Gii{p)U0,ij = I 



ipo+li- e„) - 



7oA„ 



n Pq- i<^n - m)^ - 5Tr[AjjA„](l - 63^,) + iV 

with 

A^A„ = diag (|A„p, |A„|^, O) in the color space, (32) 

where r/ is a positive infinitesimal. 

The quasiparticle energies are obtained by looking for the poles of Gii{p): 



J, , . ) V(fn(p)-M)' + |A„(?>)P for color 1, 2 

■C'n(P)= \ (33) 

V (^"(p) - A*)^ for color 3 



The quark number density Pq is also given as 



P, = -il (^Tr [(Gn(p) - Gii(p)U=o)7o] (34) 



))] (35) 



.,2/ N 1 ( . ^n{p) - H 



with 

where the first two terms in cq. (35) show the quark contributions, while the last term the anti-quark 
contribution; f^(p) is the occupation probability of the quark pairs with momentum p and represents 
difi'usencss of the momentum distribution. 

Similarly we can know the self-consistent solutions satisfy our ansatz about the mean- fields in V . From 
the above solutions we can easily obtaine that Tr[Gii(p)i75] — 0, Tr[Gii(p)7i] ocpi, Tr[Gii(p)757o] oc Pz 
and Tr[Gii(p)757i,2] oc Px,Py ■ Hence the pseudoscalar mean-field Ups, the space-component of vector 
mean-field Ul, the axial- vector mean- fields and are vanished after the integration over angles. 

C. Equation for the superconducting gap function 

Using eq. (31), the off-diagonal component of the full Green function G{p), given in eq. (8), can 
be represented in the similar way as 

G21(P) = - E ^°f2"^°iA 12^- ^n\2r2, (37) 

„ Po- (en- W^- |A„|2-|-zr/ 

where T2 is the Pauli matrix in the two-flavor space. Substituting eq. (37) into the gap equation (15) 
and using the identity X)a=i(''^a)^^2Aa = — 8/3A2, we obtain 
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7"7oB„(p)7o7" 



A„(p), 



(38) 



where the factor 2/3 is simply the Fierz coefBcient for the color and flavor degrees of freedom (Ap- 
pendix D). Furthermore multiplying both sides of eq. (38) by -B^/(fe) and taking trace with respect to 
the Dirac indices, the coupled equations for the gap functions A„ are obtained after po integration, 

^"■w = -l"7 |?5«.«('^-f)i:j-s(*.p)^ (39) 

where the function T^!^^{k,p) is defined as 

T„T„(fe,p) = TV [4,(fe)7'*7oS„(p)7o7^] = {^-n'{k)r<f>-n{p)){Mp)Y<f>n'{k)), (40) 

a decomposition of Bn{p) in terms of gamma matrices and its properties are given in Appendix B. 

Here we take the zero-range approximation in eq. (17). In terms of the polar coordinates p = 
{p, 9p,(j)p}, we can consider that the gap function A„(p) does not depend on the horizontal angle <j)p due 
to axial-symmetry around the Pz-Bxis. Thus we can explicitly perform the integration with respect to 
the angle (jjp in the gap equation (39): 



AAk,e,) = -g J -j^p sm0,^T„,„(fc,0„p,0,)^^-^ 



(41) 



with the effective coupling constant g = g/-\/Q^ + M^. As seen from the above equation, each of the gap 
functions couples with others by the function Tn'n{k,0k,p,0p) defined as 



Tn'n{k,9k,P,0p) = j ^5^,T„X(fc,p) = 



hpt 



(42) 



2\en'{k)\\en{p)\ 

where pt=psm.6p and Pz=p cos 6p and the same for kt and kz- The term proportional to Pz in eq. (42) 
will disappear after the integration over 9p. 



D. Equation for the axial-vector mean-field Ua 



Using eqs. (33) and (36), Gii(p) is recasted in the form. 



[Gii(p)]a/3, 



E 



l-vlip) 



e^f'"'A„(p)7o 



(43) 



^Po - En + ir] pq + En- iri 

Substituting the above equation into eq. (14), and integrating with respect to po, we obtain the self- 
consistent equation for Ua in the zero-range approximation: 



UA = -'-^f 

9 2^ 



(27r) 



J2[9{„-en{p)) + 2vlip)]Sn{p), 



(44) 



where the factor —2/9 stems from the Fierz coefficient of the color-singlet axial- vector channel of the 
OGE interaction (Appendix D), and Sn{p) is the expectation value of the spin operator, <Tz=—'yo7573, 
with respect to the spinor (j)n{p)- 

Ua + i-lTPp 



Snip) = Tr(7573A„(p)7o) = 4'Up){-(Tz)Mp) = 



en(p) 



(45) 



Thus Ua is related to the expectation value of az summing over the state with momentum p. An effect 
of the Cooper pairing enters into eq. (44) through the function Vn{p)- 
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E. Weak coupling approximation 



In this subsection wc consider a high-density hmit, which means the weak coupling hmit due 
to asymptotic freedom of QCD, and then disregards the anti-quark pairing and contributions from the 
negative-energy sea (the Dirac sea) in eq. (41). Actually it costs more energy to form the anti-quark 
pairing than the quark pairing for a large chemical potential. Taking the approximation, we also disregard 
the contribution from anti-quarks to calculate the quark number density in eq. (35) and the axial-vector 
mean-field eq. (44) consistently. ^ In the following calculations wc define gap functions of the quark 
pairing by subscript ± whcih corresponds to the "spin" -up (-down) of positive-energy states as A_ = Ai 
and A+ = A2. The other symbols with the subscript ± have the same meaning, e.g., 0zp = 0i,2- 

In addition, we assume that only quarks near the Fermi surface form the Cooper pairs, and thereby 
replace the gap function by an approximated form, 

A±(p)^A±(p)e(5-|e±-M|), (46) 

where (5 is a cut-off parameter around the Fermi surface. The function 6{S— \e±{p) — is also regarded as 
a form factor to regularize the integration in the gap equation [21]. The step- function form factor mimics 
the asymptotic freedom; inner particles in Fermi sea costs large kinetic energy to create the pairing 
and takes large momentum transfer which indicates that coupling of this inner-process is small. There, 
however, might be more realistic form factors for finite density QCD, which are smoother functions of 
momentum and fj. than ours, we think that they makes little change on qualitative results of the CSC 
and spin polarization. There are models with other form factors or cut-off functions [3,5]. 

Looking at the structure of the gap equation (41) with (42), one can find that the gap function is 
exactly parametrized as (Appendix C) 

e±(p) V l^P 

where R and F are some constants and represent the antisymmetric and symmetric combinations of the 
gap functions; R = /3p/m(A_ - A+) = /3p/(p(TO)(e_ A_ - e+A+) and F = A_ + A+ = l/pt(e_A_ + 
e+A_|_). Their magnitudes are determined by the coupled equations; 



I' J ^'"^^4 
2 .2 / dpd9p 2 „• /) 



Q^{p){F--R) + Q.{p){F+-R) 



r, ^ -2 / ^"P 2 • Q 



3^ y (27r)2" ^2/3p 



Tfi Tfi 

-Q+{p){F - -R) + Q-{p){F + -R) 
Pp Pp 



(48) 
(49) 



"'^^^ ^"(^)= e,(p)'k(p) ^(^-|-^(^)-^l)- 



We can obviously see that i? ^ as m ^ 0. 



^This is equivalent to the restriction of the sum over the index n{n = 1 — 4) to 1, 2, which correspond to the 
positive-energy states with different "spins" specified by the subscript =F. 
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Here we examine the polar-angle dependence of the anisotropic gap function at the Fermi surface 
A±{p^,9). The Fermi momentum p^(^) of each "spin" eigenstate is given as 



Pt = p± (0) sin 6, pz = p'^ [0) cos with 



Pi{0) = 



/i^ - + U\ cos(26l) T UA^'i^'^ cos^ + Arrfi sin^ 6-U\ sm^{20) 



1/2 



(50) 



where the subscript ± corresponds to the "spin" -up (-down) state again. Substituting the above formula 
into the gap function (47), we get 



p^{0) sin (9 



m 



m'^ + {p^{0) cose) 



-.R + F 



(51) 



Note that this form exhibits a P- wave pairing nature: it is a genuine relativistic effect by the Dirac 
spinors (Appendix B). We show a schematic view of the above gap functions in Fig. 2. 
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FIG. 2. A schematic view of polar-angle dependence of the gap functions at the Fermi surface where we set 
values of the gap parameters as i? = 0.2/i, F = //, Va = 0.3/x, (a) for m = and (b) for m = O.l/it 

As characteristic features, both the gap functions vanish at poles {9 — 0, tt) and take maximal values 
near equator {0 = 7r/2), keeping the relation, A_ > A+ ^. Suppression of A+ and enhancement of A_ 
at 6 = 7r/2 for the case of m ^ (Fig. 2b) are originated from a finite value of _R, while they vanish if 
quark is taken to be masslcss (Fig. 2a) . The anisotropic gap functions give rise to the different diffuseness 
in the momentum distribution of the two "spin" cigenstatcs. and thereby make some effects on spin 
polarization, unlike in the normal phase. The anisotropic diffuseness has two effects that it obscures the 
deformation in the momentum distribution due to their angle dependence and enlarges the difference of 
the state density between the two "spin" eigenstates through the relation A_ > A+. 



^This feature is very similar to ^P- pairing in liquid ^He [22] or nuclear matter [23]. 
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III. RESULTS AND DISCUSSIONS 



In this section we solve the coupled equations (44), (48) and (49) and investigate the effects of the 
superconducting gap on spin polarization. 

Before going to numerical calculations of Ua,R and F, each of which is coupled with others by the 
self-consistent equations, it is instructive to sec their parameter dependence by treating one of them as 
an input parameter. First we show R and F as functions of Ua in Fig. 3 where n = 400, 450 MeV and 
S = O.lAt. 




100 200 300 100 200 300 

Ua (MeV) Ua (MeV) 

FIG. 3. Parameter dependence of R and F on Ua for g = 0.13 MeV"\ m = 20 MeV and S = 0.1/u. (a) for R 
and (b) for F. Dashed (solid) lines correspond to n = 400(450) MeV. The magnitudes of R and F are calculated 
by equations (49) and (48) for given Ua- 

R starts from zero and almost linearly increases with Ua (Fig- 3a), which is understood by seeing 
that R is proportional to the difference, A_ — A+, due to finite Ua, see eq. (47) or Appendix C. Thus R 
is induced by Ua and closely coupled with it. 

As for the behavior of F, it is barely affected by Ua (slight decreasing with Ua in the numerical value) 
(Fig. 3b). As seen from the dependence on fj, the magnitude of F is almost determined by the volume of 
the phase space in the gap equation, that is, by and 6. This reflects the fact that F is related to the 
sum, A+ + A_ (Appendix C). Thus we expect that F increases with density when other parameters are 
fixed. 

From the above results we have found that F is not so much infiuenced by Ua- Next we examine the 
behavior of Ua and R when F is treated to be an input parameter. In Fig. 4 we show the parameter 
dependence of Ua and R on F, where /j. = 450 MeV and we use three values of the cut-off parameter 
d = 0.05/i, 0.1/U and 0.15/i, and add the result of Ua in the normal phase {S = 0). 
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FIG. 4. Parameter dependence of Ua and RonFfor g = 0.13 MeV"\ m = 20 MeV and n = 450 MeV. (a) for 
Ua and (b) for R. Solid line, dotted lines, dot-dashed lines and dashed lines correspond to 5 = (normal phase), 
0.05/i, O.Wfj, and 0.15/tt, respectively. The magnitudes of Ua and R are obtained by their equations for given F. 

Comparing the dependence of Ua on F (Fig. 4a) with that in the normal phase, we see a characteristic 
behavior for different values of 6: there are regions where Ua is larger than that in the normal phase 
for relatively small F, and this region seems to extend with 6. On the other hand results from the 
self-consistent calculations show that F becomes larger with S so that its value corresponds to a region 
where Ua is comparable with or slightly less than that in the normal phase, for any value of the chemical 
potential. This situation seems to be qualitatively unchanged, once the ratio of the effective coupling 
constants in the axial- vector channel G axial and the diquark channel Gdig is kept, G axial '■ Gdiq = 2/9 : 
2/3, which comes from the Fierz transformation for color and flavor (Appendix D). However, if the 
coupling constant in each channel is taken independently, our results might be changed qualitatively. 

Seeing the results for R in Fig. 4b, we find that R increases with 6 due to the growth of the phase 
space and is hardly affected by F except the region of small F where Ua varies rapidly as shown in 
Fig. 4a: it also shows that R is closely related to Ua- 

These parameter dependences also suggest that the regularization scheme for the gap equation, i.e., the 
sharp momentum cut-off function, the form factor, etc., will give rise to a qualitative change to Ua- In 
the present cut-off function, 6{6 — |e± — /i|), Ua (spin polarization) coexists with CSC, except a slight 
competition, as will be shown later. 



A. Self-consistent solutions 

We demonstrate some self-consistent solutions here. Since we have little information to determine 
the values of the parameters g and S (there may be other more reasonable form factors than the present 
cut-off function), and our purpose is to figure out qualitative properties of spin polarization in the color 
superconducting phase, we mainly set in the following calculations them as g = 0.13 MeV~^ and S = 0.1/z, 
for example, which is not so far from the couplings in NJL-like models [5,16,18]. 
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We first examine spin polarization in the absence of CSC. In Fig. 5 we show the the axial-vector 
mean-field Ua, with A± being set to be zero, as a function of baryon number density Pb{= Pq/^) relative 
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FIG. 5. Axial-vector mean-field as a function of baryon number density pu (po — 0.16 fm for g — 0.13 
MeV"^ and 5 = 0.1^. (a) for m = 14 ~ 16 MeV and (b) for m = 20 and 25 MeV. Dashed (Solid) lines are 
obtained in the normal (color superconducting) phase. 

It is seen that the axial-vector mean-field (spin polarization) appears above a critical density and 
becomes larger as baryon number density gets higher. Moreover, the results for different values of the 
quark mass show that spin polarization grows more for the larger quark mass. This is because a large 
quark mass gives rise to much difference in the Fermi seas of two opposite "spin" states, which leads to 
growth of the exchange energy in the axial- vector channel. 
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FIG. 6. R (a) or F (b) as a function of Pb/ po for m = 15 MeV (dashed lines) and m = 20 MeV (dotted lines). 
The other parameters are same in Fig. 5. Note that in figure (b) the lines almost overlap each other for the two 
quark masses. 
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Next we solve the coupled equations (44), (48) and (49). Results for Ua, R and F are shown in Fig. 5 
(solid lines) and Fig. 6, for values of the quark mass m = 14 ~ 25 MeV. It is found again, by comparing 
these cases of the quark mass, that Ua is very sensitive to the quark mass and increases with it as in the 
absence of CSC (Fig. 5). For the behavior of the gap functions, R is induced by Ua and both of F and 
R increase with pb due to the growth of the Fermi surface (Fig. 6). It is also seen that F is not sensitive 
to the quark mass. To see the bulk behavior of pairing gap as a function of baryon number density, we 
also show, in Fig. 7, their mean values with respect to the polar angle on the Fermi surface; 

The mean values (A±) begin to split with each other at a density where Ua becomes finite. This reflects 
that R is induced by Ua and then has a negative (positive) contribution to A+ (A_). 




2 4 6 8 10 12 14 

Pe/Po 

FIG. 7. Mean values of the gap functions with respect to the solid angle at the Fermi surface, (A±), plotted 
as a function of ps/po for m = 15 MeV (solid and dotted lines) and 20 MeV (dot-dashed and dashed lines). 

Here we would like to comment on the magnitude of (A±). These should be compared with the 
usual uniform gap function, and may look very large values of O(GeV) in our case. However these values 
would be largely reduced by taking a smooth form factor which models asymptotic freedom of QCD [5] ; 
it further reduces the integral value in the gap equation, compared with our sharp cut-off function. 

In Fig. 8 wc show the expectation value of the spin operator per quark, (cr^/A^q), as a function of 
Pb/po with and without the superconducting gap. The critical density becomes lower as the quark mass 
increases, and the peak positions of [uz/Nq) are located at relatively lower densities in each quark mass. 
The magnitude of (a^/Nq) is to be compared with 1 for a free quark, because \tljl(7zips/'4'l'4's \ = 1 at the 
rest frame for the free spinor ips- We arrange the results of three quark masses m = 14 ~ 16 MeV by 1 
MeV in Fig. 8a to show a high sensitivity of spin polarization to the quark mass, which implies that the 
exchange energy from the attractive axial-vector interaction is strongly enhanced by the quark mass to 
produce the large axial-vector mean-field. 

The exchange energy is also enhanced by larger chemical potential and the resulting axial-vector 
mean-field increases with it (see Fig. 5). But the spin expectation value per quark, which is relative to the 
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axial-vector mean-field per quark (a UA/Nq), has an upper limit since the increase of Nq is far superior 
to that of Ua for larger chemical potential, which gives rise to the peak positions in Fig. 8. 
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FIG. 8. Spin expectation value per quaxk as a function of Pb/po- (a) for m = 14 ~ 16 MeV and (b) m = 20 
and 25 MeV. Dashed (Solid) lines show results in the normal (color superconducting) phase. 

The quark mass is very important in relation to the breaking of chiral symmetry in QCD. Models 
incorporating chiral dynamics have indicated that the dynamical mass becomes smaller as chiral symmetry 
is restored at a high density, while the current quark mass is small and explicitly breaks it [24]. In our 
model, on the other hand, we treat the quark mass m as a variable parameter so that we may simulate 
a change of the dynamical mass. In order to further examine the effect of the quark mass on spin 
polarization, we show the mass dependence at densities ps = 5po, pB = lOpo and ps = 15po for the 
cases with and without the superconducting gap in Fig. 9. 
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FIG. 9. Spin expectation value per quark as a function of the quark mass m for fixed baryon number density 
ps/po = 5, 10 and 15. 
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Spin polarization increases with the quark mass in all the three densities. In the figure we exhibit only 
the results for a narrow region of the mass parameter (m = 13 20 MeV), while as for larger masses of 
O(lOOMeV) (order of the strange quark mass) spin polarization monotonically increases without singular 
oscillations. Critical values of the quark mass at which spin polarization disappears become smaller as 
density increases in both cases. 

In relation of Ua to m we can derive an exact result in the massless limit, m — > 0. In the normal 
phase where A = 0, eq. (44) becomes 

^^ = -q9 l^J (^30(M-e„(p)) (53) 

with 



e±(p) = ^{\p.\±UAf +pI (54) 
The right-hand side of the above equation can be analytically integrated to give 



Ua = --^9 jT^z^-i / / dptpt 

Jo Jo 



r (27r)3 j, ^^"Vb. + W+Pt^ 

/o / (ip^ / dptpt- 



9 (27r)3 Jo Jo ViPz-UA)^+Pt 

= (55) 

Here we have assumed that > Ua- In the massless limit, the Fermi sea is described by two complete 
spheres in the momentum space with radii fj,, whose centers are located at {pt,Pz) = (0, ±Ua) (see Fig. Ic). 
The momentum distribution for quarks in the "spin" -down state occupies these two spheres, while the 
"spin" -up state does their overlap region (shaded). In the above integration for the "spin" -down quarks, 
the expectation value of the spin operator given by quarks with < Pz < Ua is canceled with that of 
quarks with Ua < Pz < "^Ua- The remaining contribution from the region, 2Ua < Pz < tJ- + Ua, cancels 
with that by the "spin" -up quarks. Thus we can see that spin polarization disappears as m — > in the 
absence of CSC. 

This analytical result that Ua — > as to ^ can be also understood as follows. The eigenstates of 
non-interacting massless fcrmions arc classified by the definite lic^licity states: the positive energy state 
is right-handed (left-handed) with positive (negative) helicity, while the negative energy state those with 
negative (positive) helicity. This property is not spoiled by introducing the axial- vector mean-field, when 
we extend the meaning of helicity; the Dirac equations for the "left-" and "right-handed" positive-energy 
fermion fields ipL,R are now given as 

{po+p-cr + UA(T3)tpL = (56) 
{pQ-p-cr + UA(T3)tpR = 0, (57) 



which give the eigenvalues, pa = yjpj + {pz + Ua)^{= £l{p)) for ipL and po = + (P^ ~ Ua)^{^ ^r{p)) 
for iliR. xpL O^r) is the eigenstate of generalized helicity h — ^1 projected onto the shifted momentum 
p' = {Px,Py,Pz ± Ua}- If /X ^ they form the spherical Fermi seas, sec Fig. Ic. Here it would be 
interesting to compare these eigenvalues with the limit form of e±{p) in eq. (22), 
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Then we can see the relations: e±(p) = eL{p)0{±Pz) + £r{p)0{TPz), which clearly show that the two 
Fermi seas of the eigenspinors (23) give the same Fermi seas of ipL.n. for a given chemical potential fi. Thus 
we can take an alternative view of the Fermi seas in terms of the definite hclicity states by rearranging 
the cigenspinor (23) properly in the massless limit. In each Fermi sea for iJjl.r the particle number with 
the definite h becomes the same, and thereby the total spin-expectation value becomes vanished. 

In the color superconducting phase, on the other hand, the situation is different because the momen- 
tum distribution becomes diffused due to the creation of the Cooper pairs near the Fermi surface. For 
m — » and A± ^ 0, eq. (44) becomes 

2, 



(27r) 



en(p) 



(59) 



where v^{p) indicates the diffused part of the momentum distribution, defined in eq. (36). Here we should 

note that the gap functions A± are still non zero even at m = 0. The diffused part, however, give no 
contribution to the spin polarization from the viewpoint of the helicity eigenstates. The gap function in 
the massless limit becomes 



pt 



(60) 



V(b.|± W+p^ ' 

see Fig. 2a. The diffuseness from the above gap function has an equal contribution to the two complete 
Fermi spheres of chirality. Thus the total spin, which is obtained by summing up these momentum 
distributions, should be zero in the massless limit even if the CSC is taken into account. 

To summarize we show a phase diagram for the quark mass and baryon number density in Fig. 10 
where we add the result of g = 0.26MeV~^ to see the dependence on the coupling constant. 
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FIG. 10. Phase diagram in ps/po — m plane for the effective coupling constant g = 0.13 MeV"'^ and 0.26 
MeV"'^. At regions above the lines spin polarization arises. This result is obtained with the gap function, thus 
all the region in the phase diagram shows the color superconducting phase. 



18 



The lines indicate the critical mass at a fixed density (at regions above the lines spin polarization 
arises). We can confirm that the critical mass becomes smaller with the increase of the density, and spin 
polarization occurs at moderate densities (ps = 3 4po) if the coupling is strong enough even though 
quark mass is taken to be smaller as a simulation for change of dynamical mass (restoration of chiral 
symmetry). 

Here we would like to understand how the gap function affects spin polarization and brings about 

a slight reduction of it. In the spin-polarized phase, the momentum distribution is deformed from the 
simple spherical shape. As mentioned in Ref. [13], the deformation is induced by finite Ua and feeds 
back to Ua in a self-consistent manner. In the color superconducting phase, diffuseness caused by the 
Cooper pairing in the momentum distribution depends on the polar angle and then has an influence on 
the deformation. As can be expected from the polar-angle dependence of the gap function, diffuseness 
tends to obscure the deformation. 

From the consideration of the spin expectation values by spinors (23) near the Fermi surfaces; 

0^(-CT^)0± = ^- W , (61) 

where (pzp = 4>i.2 for two "spins". The difference of the spin expectation value between two "spin" states 
(j)± is largely affected by high-p^ regions or regions near both poles {0 = 0, tt). Thus the large deformation 
along the z-axis seems to enhance spin polarization. 

In order to specify to what extent the Fermi sea is deformed, we calculate the quadrupole deformation 
of the momentum distribution defined by 

Q2 ^ 3(p^)/(p2) _ 1. (62) 

In Fig. 11, we show Q2 as a function of Ua at /i = 450 MeV in the normal phase and in the color 
superconducting phase in which the gap functions are given by their equations for fixed Ua- 




0.0 0.1 0.2 0.3 0.4 0.5 

Ua/a^ 

FIG. 11. Quadrupole deformation of the momentum distribution (62) as a function of Ua- Parameters are 
fixed as g = 0.13 McV~\ m = 20 MeV and n = 450 MeV. Dashed (solid) line is given in the normal (color 
superconducting) phase. 
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Prom this result of Q2 deformation, we can see that the diffused part near the Fermi surface obscures 
the deformation then gives an opposite effect against Q2, and thus reduces spin polarization. 

Nevertheless the gap function has another effect on spin polarization. It is to be noted that the 
qualitative relation, A_ > A+, is always retained as seen from Fig. 2 and then has a effect to enlarge 
the difference of the state density between the two "spin" states. This effect is expected to enhance 
spin polarization since the difference of the spin expectation value by each spinor (61) near the equator 
(6 = 7r/2), so that pz ~ 0, seems to depend only on the difference of the state density. To see it in both 
the normal and color superconducting phases, we define that Nup {Ndown) is the state density of the 
"spin"-up (-down) state and show their difference by dN = Ndown — Nup, only in the first two colors, as 
a function of Ua in Fig. 12 at = 450 MeV. The result indicates that the gap functions slightly enhance 
dN than normal phase. 
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0.0 0.1 0.2 0.3 0.4 0.5 

FIG. 12. Difference of the state densities in the two "spin" states plotted as a function of Ua- Legends are the 
same as in Fig. 11. 

From the above discussions spin polarization is significantly influenced by both the deformation and 
the state density in each "spin" state. As a result of the self-consistent calculation in the color super- 
conducting phase, the reduction effect on the deformation is slightly superior to the enhancement effect 
from the difference of the state densities, and the pairing effect finally reduces spin polarization than in 
the normal phase. It, however, should be noted that this qualitative conclusion about whether CSC en- 
hances spin polarization than normal phase or not is very delicate and may be changed depending on the 
regularization scheme, as already mentioned. Moreover other types of pairing which are not considered 
here, e.g. pairing of the "spin" -up and -down states, may gives rise to qualitatively different results, 
while it is very difficult to see which type of pairing is energetically favored. 

Finally we would like to comment on the coupling of the spin polarized quark matter with the 
external magnetic field; quark fields couple with the magnetic field through its anomalous magnetic 
moment. The magnetic interaction is described by the Gordon identity for the gauge coupling term: 
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gLe* /2m{'ipanvtp)F^^'^ where ql is a form factor and e* an effective charge. ^ In quark matter a magnetic 
moment is given as the expectation value {aij) with respect to the ground state. In our model only {(Tu) 
is nonzero, and the magnetic moment per quark is given as 

= (ai2/iV,) = -J2 I 7?T3 + ^(^^ " en)]^n(p)ai20„(p). (63) 

Note that the expectation value of <ti2 by the spinor does not depend on Ua; ^±(p)o'i20±(p) = T^/Pp, 
so that Mz reflects only the asymmetry in the momentum distribution due to the axial- vector mean-field. 




FIG. 13. Induced magnetic moment per quark (63) as a function of ps/po- Parameters and legends are the 
same as in Fig. 5. 



In Fig. 13, Mz is given as a function of baryon number density. This indicates that resulting ground 
state also holds ferromagnetism (spontaneous magnetization). 

IV. SUMMARY AND CONCLUDING REMARKS 

In this paper we have examined spin polarization in quark matter in the color superconducting 
phase. We have introduced the axial- vector self-energy and the quark pair field (the gap function) , whose 
forms are derived from the one-gluon-exchange interaction by way of the Fierz transformation under 
the zero-range approximation. Within the relativistic Hartree-Fock framework we have evaluated their 
magnitudes in a self-consistent manner by way of the coupled Schwinger-Dyson equations. 



*Here we needn't consider the orbital angular momentum for uniform matter. But if a superconductor is of the 
'second' type in which London's penetration length is larger than the coherence length, a vortex lattice may be 
formed in response to the external field and then total magnetization is to undergo a qualitative change due to 
circulation of supercurrent [25] . 
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As a result of numerical calculations spontaneous spin polarization occurs at a high density for a finite 
quark mass in the absence of CSC, while it never appears for massless quarks as an analytical result. In 
the spin-polarized phase the single-particle energies corresponding to spin degrees of freedom, which are 
degenerate in the non-interacting system, are split by the exchange energy in the axial-vector channel. 
Each Fermi sea of the single-particle energy deforms in a different way, which causes an asymmetry in 
the two Fermi seas and then induces the axial-vector mean-field in a self- consistent manner. In the 
superconducting phase, however, spin polarization is slightly reduced by the pairing effect; it is caused by 
competition between reduction of the deformation and enhancement of the difference in the phase spaces 
of opposite "spin" states due to the anisotropic diffuseness in the momentum distribution. 

In connection of the deformation with superconductivity it has recently been reported [26] that in the 
superconducting asymmetric nuclear matter the Fermi sea may undergo a deformation even in the spin- 
saturated system due to the difference of the Fermi surface between neutrons and protons; the momentum 
distributions of neutrons and protons may deform respectively to enlarge the overlapped region in the 
phase space, which effectively contributes to the np- pairing. They have shown the possibility of the 
deformation in a variational way; the Fermi sea of the majority of nucleons deforms in a prolate shape, 
while the minority in an oblate shape. Thus the deformation property of the Fermi seas looks very similar 
to our case. Nevertheless, note that our deformation is produced by the relativistic effect. Anyway it 
would be interesting to look further into the common feature. 

It is to be noted that if the effective coupling constant is strong enough to lower the critical quark 
mass, spin polarization (magnetization) has potential to appear at rather moderate densities such as in 
the core of neutron stars, even though CSC weakly works against it. 

From the above observations it is suggested that spin polarization does not compete with CSC but can 
coexist with it, unlike in ordinary superconductors of the electron system with the s-wave and spin-singlet 
pairing. This reflects the fact that internal degrees of freedom of the quark field, e.g.the color, flavor and 
Dirac indices, have rich structures to satisiy the antisymmetric constraint on the quark-pair field. 

The possibility of the coexistent phase might also give a clue as for the origin of the superstrong 
magnetic field observed in magnetars. We roughly estimate the expected magnetic field when magnetars 
are assumed to be quark stars. The maximum dipole magnetic field at the star surface reads 

Smax = —^gn,{{M,)/Ng), (64) 

with fj,g and Uq being the quark magnetic moment and the quark number density, respectively;e.g., for 
{M,)/Ng - (10-3) and ~ ©(Ifm'^), we find Bmax ~ O(IO^^G), which is comparable to that observed 
in magnetars [8,9]. 

In the present paper we have not taken into account chiral symmetry, which is one of the basic concepts 
in QCD. If chiral symmetry is restored at finite baryon number density, the quark mass becomes drastically 
smaller as density increases. In order to simulate it we have examined the quark mass dependence on 
spin polarization. In the future work we would like to consider an effect of the dynamical mass on the 
axial- vector self-energy. 
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APPENDIX A: STRUCTURE OF SPINOR UNDER THE AXIAL- VECTOR MEAN-FIELD Ua 



In this Appendix we rewrite the spinor (23) in terms of the free quark one and the remainder charac- 
terized by Ua- We employ the free spinor Us{p) in which the two-component Pauli spinors are given as 
eigenvectors of the spin matrix a^: 



u±{p) 



with C+ = I and ^_ = 

/ I 1 



(Al) 



where eo = ■\/p'^ + and {cr} are the Pauli spin matrices. 

The spinor 4'+{p) = 02 (f>) for the "spin" -up state is decomposed as follows : 
/3p(Ae - 2Ua) + m{5e + 2/3^) 



-Mp) 



ofl / ] ,'e+-Pp-UA , . Pp + m . A , Remi([/A) 

2ppVeo + m ( — - — u+{p) ^ W-(p) ) + 



Px+iPy Pz "7 ' PziPx+iPyV 

where = e_ — e+, Ae = e_ -|- e+ and 

/ p,{e+ -Pp- Ua) [PpiAe - 2Ua - 2eo) + m6e] 

-iPx + iPy){Pp + m) [f3p{Ae - 2Ua - 2eo) + mSe] 
5e[{e+ -/3p- Ua)pI + 2/3^p{Pp + m)] - 2/3p(/3p + m) [(e_ - UA){e+ - Ua) - eg] 
V Pz {Px + iPy ) (/3p + m) (5e 



(A2) 



Remi(?7A) = 



(A3) 



Note that the term Remi({7A) vanishes in the limit, Ua — > 0. Thus one can fined that <p+{p) is a mixture 
of the free spinors even when Ua = 0- 

A decomposition of the spinor 0-(p) = (f)i{p) for the "spin" -down state can also be done in the similar 
way: 

PpjAe - 2Ua) + m{6e + 2pp) 



PtPz-M'- 



-ip) = 



2/3p\/eo + m 



e-+Pp-UA , . (3p-m_ ^ Rem2(C/A) 



Px + iPy 



-u+{p) + 



-u-ip) 



Pz{Px +iPy)' 



(A4) 



where 



Rem2(f/A) = 



/ p,{e- +Pp- Ua) [/3p(Ae - 2Ua - 2eo) + mSe] 

(p^ + ipy){Pp - m) [l3p{Ae - 2Ua - 2eo) + mSe] 
Se[{e- +(3p- Ua)pI - 2Pl{Pp - m)] + 2(3p{Pp - m) [(e_ - [/^)(e+ - Ua) - eg] 
V -Pz{Px + iPy){l3p - m)6e 



(A5) 
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APPENDIX B: DECOMPOSITION OF S„(P) IN TERMS OF THE DIRAC GAMMA 

MATRICES 



The operator Bn{p) in eq. (27) consists of some gamma matrices; it is a linear combination of 1, 7, 
757) CToi and cro2, in the diquark field ipcBntp = ^j'^CBnip. The last two matrices give the tensor diquark 
fields, while these terms have no influence on the gap equation (38) due to axial symmetry of the Fermi 
seas around the Pz axis: the integration of Bn{p) with respect to the azimuthal angle (/>p in eq. (38) gives 

Bn{p) ^ r ^Br^ip) = f [(-l)>z73 + (-l)"ml + /?p7573] • (Bl) 
Jo 27r ^en[p)\Pp 

Thus tensor terms disappear because they are proportional to exp (i^p) in Bn(p). 

The first term in the right hand side has also no contribution after symmetric integration with respect 
to pz- The remainders, {1,7573}, imply the pseudo-scalar (J^ = 0^) and vector (J^ = 1^) diquark 
pairings in terms of the notation in ref. [2]. Please note that the CSC gap (Bl) results in a linear 
combination of different angular momentum pairs 0~ and 1~ because of the lack of rotation symmetry. 

Since the diquark fields ^^(7(1,7573)^ contain the off-diagonal matrices which connect the lower 
component with the upper one of the Dirac spinors, these pairings vanish in the non-relativistic limit 
or in the limit m — > 00. Hence -Bn(p) resembles P- wave pairing as is seen in eq. (50), although it has 
no correspondence in the non-relativistic limit: the gap function for (Bl) has the nodes (vanishing at 
^ = 0, 7r) due to the factor pt, which is similar to ^P-pairing in the liquid ^He - A phase, but these nodes 
are entirely attributed to the genuine relativistic effect. This property survives even in the limit, Ua 0. 

From eq. (Bl) we can also obtain the relation appearing in eq. (40); 

7/*7o^„(p)7o7'' = 2S„(p) +2m{m+ (-l)"/?„} . (B2) 



APPENDIX C: PARAMETERIZATION OF THE GAP FUNCTION 



In this Appendix we derive the parameterization (41). The gap equation (47) is expanded as 

Pt (, 2rn2 \ A^(p) jp^_ / 2m^ \ A_(p) ' 
e+(p) V hdp J ^E+ip) ^ e.{p) V(3k(3p ) 2E_{p)_ 



^^('^ = 3^ i (2^2;^ 



Introducing A± (k) through the equation, 



(Cl) 



(C2) 



we obtain the "gap" equation for A±{k), 



A m = 2.2 f^pi 2m^ / A.jp) A+{p) \ A_(p) A+jp) 

' 3^ J (27r)3 4 [e.{pyE_{p) e+{pYE+{p) ) ^ e.{pYE_{p) ^ e+{pYE+{p) _ 



Then we find the following properties. 



A_{k) + A+{k) = F and A_{k) - A+{k) = R 



/3k' 



(03) 



(C4) 
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where F (i?) is a constant which characterize the symmetric (asymmetric) combinations of the gap 
functions A±. Thus we can further parameterize As(fc) as 

777 

A±{k)=T^R + F. (C5) 

Pk 

Substituting the above formula into eq. (C3), one can obtain the coupled equations for F and R, eqs. (48) 
and (49). 



APPENDIX D: FIERZ TRANSFORMATION 



Wc present the Fock exchange energy term by the OGE interaction by the use of the Fierz transfor- 
mation. The Green function with vertices in the right-hand side of eq. (14) can be expanded as 

^ (r„iGn(p)r„),. = Y.^Ta)iAmi-mj'){^a)j'j = 5^a6(r6)i,Tr(Gnr6) (Dl) 

a a ab 

with Ta = ^u.® ^flavor <S) Kolor and {Ta)ii'{Ta)j'j = ^ Cab{Yb)ij{^b)3'i' , (D2) 

6 

where {Cab} are coefficients of a Fierz transformation (D2) for the Dirac matrices, the identity matrix in 
the flavor space and the Gell-Mann matrices in the color space; 



{lti)iv{l'')j'j = - - ^(757/x)y(757^)i'z' + (*75)*j (ns)^' 



Sii'Sj'j — 



1 



■SijSj'i' + {Ta)ij{Ta)j 



1 



i^c)ii'iK)j'] — "j^^i^c ^ l)%<5j'i' - ^(^c)ii(Ac)j'i' 



(D3) 
(D4) 
(D5) 



It is to be noted that there appears no tensor term in eq. (D3) due to chiral symmetry in QCD. Thus, 
e.g., the coefficient for the color-singlet axial- vector self-energy reads —4/9 for Nf ~ 2 and Nc = 3. 

We also present a Fierz transformation for diquark fields. The right hand side of eq. (15) can be 
expanded, in the similar way for Gn; 

^(f„G2i(p)r„),,- = J2iCTlc-%,{Mp)Mp)rK^ah'i = E('^)^fe(ra)i'fc(Vi(-p)i'V;(p),')(ra)i',- 

a a a 

= ^fabTT {G2iip)C-'Tl C-') (CTi'C^j (D6) 

ab 

with (r„)i,fe(r„),v,. = ^/„(,(r6C*)i.,v(CT6),fe, (D7) 

b 

where {fab} are coefficients of a Fierz transformation (D7) and are explicitly given as 



Si'kSj'j 



5i'jiSjk + {Ta)i'j'{Ta)jk 



(Ac)i'fe(A, 



(D9) 
(DIO) 



where {Af^"^"*} are symmetric (antisymmetric) matrices of {Ac}. The present gap function, A(p) = 
J2n -Bn(p)A„(p) which is a linear combination of the gamma matrices, can be obtained by taking pro- 
jection on i?„(p). 
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